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We derive Ginzburg Landau (GL) expansion in disordered attractive Hubbard model within the 
combined Nozieres - Schmitt-Rink and DMFT+E approximation. Restricting ourselves to the case 
of homogeneous expansion, we analyze disorder dependence of GL expansion coefficients on disorder 
for the wide range of attractive potentials U, from weak BCS coupling region to the strong coupling 
limit, where superconductivity is described by Bose - Einstein condensation (BEC) of preformed 
Cooper pairs. We show, that for the case of semi - elliptic “bare” density of states of conduction 
band, disorder influence on GL coefficients A and B before quadratic and fourth - order terms of 
the order parameter, as well as on the specific heat discontinuity at superconducting transition, is of 
universal nature at any strength of attractive interaction and is related only to the general widening 
of the conduction band by disorder. In general, disorder growth increases the values of coefficients A 
and B, leading either to the suppression of specific heat discontinuity (in the weak coupling limit), 
or to its significant growth (in the strong coupling region). However, this behavior actually confirms 
the validity of the generalized Anderson theorem, as disorder dependence of superconducting critical 
temperature T c , is also controlled only by disorder widening of conduction band (density of states). 


I. INTRODUCTION 


The problem of superconductivity in BCS — BEC 
crossover region (and up to the strong coupling limit) has 
a long history, starting with early works by Leggett and 
Nozieres and Schmitt-Rink E,a. Probably the simplest 
model to study this crossover is Hubbard model with at¬ 
tractive interaction. The most successive approach to 
the studies of Hubbard model (both repulsive and at¬ 
tractive) is the dynamical mean field theory (DMFT) 
Attractive Hubbard model was already studied 
within DMFT in a number of papers (M3. However, 
up to now there are only few works, where disorder ef¬ 
fects were taken into account, either in normal or super¬ 
conducting phases of this model. Qualitative analysis of 
disorder effects upon critical temperature T c in BCS — 
BEC crossover region was presented in Ref. PI, which 
claimed the validity of Anderson theorem in this region 
for the case of s-wave pairing. Diagrammatic analysis of 
disorder effects on T c and the properties of the normal 
state in crossover region was recently presented in Ref. 


We have developed the generalized DMFT+E ap¬ 
proach to Hubbard model (l3lT16j . which is quite conve¬ 
nient for the account of different “external” interactions, 
e.g. such as disorder scattering [17, JJ1- This approach 
is also well suited to the studies of two-particle prop¬ 
erties, such as dynamic (optical) conductivity pi, |19j . 
In a recent paper 0 we used this approach to ana¬ 
lyze the single-particle properties of the normal phase 
and optical conductivity in attractive Hubbard model. 
Further on the DMFT+E approximation was combined 
with Nozieres - Schmitt-Rink approach to study the in¬ 
fluence of disorder on superconducting critical temper¬ 
ature T c in BCS BEC crossover and strong coupling 


region [20l . l 2l| , demonstrating the validity of the gener¬ 
alized Anderson theorem. Disorder effects upon T c are 
essentially due only the general widening of the conduc¬ 
tion band by random scattering. This was demonstrated 
exactly (for the whole range of attractive interactions) 
in the case of semi - elliptic density of states of conduc¬ 
tion band (three-dimensional case) at any disorder level 
and becomes also valid in the case of flat band (two- 
dimensional case) in the limit of strong enough disorder. 

Ginzburg - Landau (GL) expansion in the region of 
BCS BEC crossover was derived in a number of previous 
papers [22|-|24j], however no effects of disorder scattering 
on GL - expansion coefficients was considered. Here we 
derive the microscopic coefficients of (homogeneous) GL 
- expansion for the attractive Hubbard model and study 
disorder effects on these coefficients including the BCS 

BEC and strong coupling regions, as well as upon the 
specific heat discontinuity at superconducting transition, 
demonstrating certain universality of disorder behavior 
of these characteristics.. 


II. DISORDERED HUBBARD MODEL IN 
DMFT+E APPROACH 

We consider the disordered attractive Hubbard model 
with Hamiltonian: 

H = -i X a «r a U + X einia n *t n 4> (!) 

(ij)cr icr i 


where t > 0 is transfer integral between the nearest neigh¬ 
bors and U is onsite Hubbard attraction , = a\ a a i(T is 

electron number operator at site i, a,i a (aj^) is annihila¬ 
tion (creation) operator of an electron with spin a. Local 
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energy levels are assumed to be independent random 
variables on different lattice sites. We assume the Gaus¬ 
sian distribution of ti at each site for the validity of the 
standard “impurity” scattering diagram technique |25j : 

= (2) 

Here A is the measure of disorder scatte ring. 

The generalized DMFT+E approach |l3fil6l | supplies 
the standard DMFT [3h5j with an additional “external” 
self-energy (in general case momentum dependent) due 
to any interaction outside the DMFT, which provides 
an effective method to calculate both single and two - 
particle properties 0 , 0 . The additive form of the to¬ 
tal self-energy conserves the structure of self - consistent 
equations of DMFT jSH'SS- The “external” self-energy is 
recalculated at each step of the standard DMFT itera¬ 
tion scheme, using some approximations, corresponding 
to the form of an additional interaction, while the local 
Green’s function (central for DMFT) is also “dressed” by 
additional interaction. 

For disordered Hubbard model we take the “external” 
self-energy entering DMFT+E cycle in the simplest form 
of self - consistent Born approximation, neglecting the 
“crossing” diagrams due to disorder scattering: 

E( £ ) = A 2 ^G(+p), (3) 

p 

where G( £ , p) is the complete single - particle Green’s 
function. 

To solve the effective Anderson impurity model of 
DMFT here we used the effective alg orithm of numer¬ 
ical renormalization group (NRG) |2q |. 

In the following, we consider the model of the “bare” 
conduction band with semi - elliptic density of states (per 
unit cell and spin projection): 

N o(e) = ^V D2 -e* (4) 

where D determines the half - width of conduction band. 
This is a good approximation for the three - dimensional 
case. 

In Ref. [~21] we have shown analytically, that in 
DMFT+E approach, within these approximations, all 
the disorder influence upon single - particle properties 
is reduced to the simple effect of band widening by dis¬ 
order scattering, so that D —> D e ff , where D e ff is the 
effective band half - width in the presence of disorder (in 
the absence of correlations, i.e. for U = 0): 

Deff = D^J 1 + 4^. (5) 

and conduction band density of states (in the absence of 
U ) “dressed” by disorder is given by: 

«o( £ ) = (6) 

71U eff 


conserving its semi - elliptic form. 

For other models of the “bare” conduction band den¬ 
sity of states, besides band widening, disorder scattering 
changes the form of the density of states, so that the com¬ 
plete universality of disorder influence of single - particle 
properties, strictly speaking, is absent. However, in the 
limit of strong enough disorder the “bare” band density 
effectively becomes elliptic for any reasonable model, so 
that universality is restored (2l|. 

All calculation below were performed for the quarter - 
filled band (n=0.5 per lattice site). 


III. GINZBURG LANDAU EXPANSION 

The critical temperature of superconducting transition 
T c in attractive Hubbard model was analyzed using direct 
DMFT calculations a number of papers @,0i. In Ref. 
01 we have determined T c from instability condition of 
the normal phase (instability of DMFT iteration proce¬ 
dure). The results obtained in this way were in good 
agreement with the results of Refs, [a, 13, 01- Addition¬ 
ally, in Ref. 0 we calculated T c using the approximate 
Nozieres - Schmitt-Rink approach in combination with 
DMFT (used to calculate the chemical potential of the 
system), demonstrating that being much less time con¬ 
suming, it provides semi - quantitative description T c 
behavior in BCS - BEG crossover region, in good agree¬ 
ment with direct DMFT calculations. In Refs. [2(j, [21) 
the combined Nozieres - Schmitt-Rink approach was used 
to study the detailed dependence of T c on disorder. Be¬ 
low we shall use this combined approach to derive GL 
- expansion including the disorder dependence of GL - 
expansion coefficients. 

We shall write GL - expansion for the difference of 
free energies of superconducting and normal phases in 
the standard form: 

F s -F n = A|A q | 2 + q 2 C |A q | 2 + ||A q | 4 , (7) 

where A q is the spatial Fourier component of the am¬ 
plitude of superconducting order parameter. Microscop¬ 
ically, this expansion is determined by diagrams of the 
loop - expansion for the free energy of an electron in the 
“external field” of (static) superconducting order parame¬ 
ter fluctuations with small wave vector q , shown in Fig'll] 
(where fluctuations are represented by dashed lines) |25| . 
Below we limit ourselves to the case of homogeneous ex¬ 
pansion with q = 0 and calculations of its coefficients A 
and B , leaving the (much more complicated) analysis of 
the general inhomogeneous case of finite q and calcula¬ 
tions of coefficient C in (0 for the future work. 

Within Nozieres - Schmitt-Rink approach 0 we use 
the weak coupling approximation to calculate loop - di¬ 
agrams with two and four Cooper vertices shown in Fig. 
IH dropping all corrections due to Hubbard U, while in- 
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Using the symmetry e(p) = e{— p) and G(e n ,— p) = 
G(e„, p), we obtain from the Ward identity (fTTH) : 


y ] ^pp' ( £ n) — 

pp' 


E P G(e n , p) - E P G(-e n , p) 


2 is 


n 


(ii) 


so that for Cooper susceptibility © we get: 


Xo(q = 0 ,T) = 

_ T Ep G(£n, p) - Ep G(-£ n , p) 

n 2l£n 

T Sp Epg(£n,p) 


( 12 ) 


FIG. 1: Diagrammatic representation of Ginzburg Landau Performing now the standard summation over Matsubara 
expansion. frequencies (2fj, we obtain: 


eluding “dressing” by disorder scattering 1 . However, the 
chemical potential, which essentially depends on the cou¬ 
pling strength U and determines the condition of BEC 
in the strong coupling region, is calculated via the full 
DMFT+E procedure. 

Coefficient A before the square of the order parameter 
in GL - expansion is given by diagrams of Fig. [TJ a) with 
g = 0@: 


A(T) = X o(q = 0, T) - X o (« = 0, T c ), (8) 

where 

X o(5 = 0,T) = -T^5> pp ,( e „) (9) 

n pp' 

is the two - particle loop in Cooper channel “dressed” 
only by disorder scattering, while 'bpp, (e n ) is disorder 
averaged two - particle Green’s function in Cooper chan¬ 
nel (e n = 7rT(2n + 1) is corresponding Matsubara fre¬ 
quency). Subtraction of the second diagram in Fig. [fja), 
i.e. that of Xo(q = 0, T c ) in (|5|). guarantees the validity 
of A(T = T c ) = 0, which is necessarily so in any kind of 
Landau expansion (25j . 

To obtain X) PP ' ^pp'( £ n) we use an exact Ward iden¬ 
tity, derived in Ref. Qj)|: 


Xo {q = 0 ,T) = 


1 

47 Tl 


7-00 £ 22 



, N(s) , £ 

de - -th —, 

2e 2 T’ 


(13) 


where N(e) is the “bare” (U = 0) density of states, 
“dressed” by disorder scattering, which in the case of semi 
- elliptic band takes the form ©. In Eq. (fT3l) the origin 
of e is at the chemical potential. Replacing s —> e — /r 
to move the origin of energy to the center of conduction 
band, we finally write: 


Xo {q = 0, T) 



de 


N(s) 

2(e - n) 


th 


£ — H 

2T ’ 


(14) 


Cooper instability of the normal phase, determining su¬ 
perconducting transition temperature T c , is written as: 


1 = —U Xo(q = 0, T c ) (15) 


Then, to determine the critical temperature we obtain 
the following equation: 




deNo{e) 


th^_ 

£ — fl ’ 


(16) 


G(£ n ,P) - G(-£ n , -p) = 

= -^^pp'^")^ 1 ^’?') - GQ 1 (—e n , —p')) ? (10) 

P' 

Here G(e n , p) is disorder averaged (but not “dressed” by 
Hubbard interaction!) single - particle Green’s function. 


1 In the absence of disorder this approach just coincides with that 
used in Refs. (234231 , using Hubbard - Stratonovich transforma¬ 
tion in the functional integral over fluctuations of superconduct¬ 
ing order parameter. 


Using m to determine Xo(q = 0, T c ) and (ThO) for X o(q = 
0,T), we obtain the coefficient A (jSJ): 

The chemical potential for different values of U and A 
is to be determined here from direct DMFT+E calcu¬ 
lations, i.e. from the standard equation for the total 
number of electrons (band filling), defined by Green’s 
function obtained in DMFT+E approximation. This al¬ 
lows us to find both T c and GL - expansion coefficients 
in the wide range of parameters of the model, including 
















4 


the BCS - BEC crossover region and the limit of strong 
coupling, for different disorder levels. Actually, this is the 
essence of Nozieres - Schmitt-Rink approximation — in 
the weak coupling region transition temperature is con¬ 
trolled by the equation for Cooper instability, while in 
the strong coupling limit it is defined as the temperature 
of Bose condensation, which is controlled by the equation 
for chemical potential. The joint solution of Eqs. m 
and m with DMFT+E equation for chemical potential 
provides the correct interpolation for T c and GL - coef¬ 
ficient A from weak coupling region via the BCS - BEC 
crossover towards the strong coupling. 

For T —> T c the coefficient A(T ) is written as: 

A(T) = a(T — T c ). (18) 


where in case of temperature independent chemical po¬ 
tential 


a = 


1 

4 T* 



<IeNq(e) 


ch*& 


(19) 


In BCS approximation with conduction band of infinite 
width with constant density of states No(0) we obtain 

from (fT9l) the standard result a = [25.]. However, in 

BCS - BEC crossover region temperature dependence of 
/r is essential and we have to use the general expression 
m in conjunction with equation for /r to calculate a. 
At the same time, from Eq. E3) it is clear that disorder 
scattering influences a only through the changes of the 
density of states Nq(e) and chemical potential //, which 
is the typical single - particle property. Thus, in the case 
of semi - elliptic “bare” conduction band the dependence 
of a on disorder is due only to the band widening by dis¬ 
order replacing D —> D e ff. Thus, in the presence of dis¬ 
order we expect the universal dependence of a(2D e ff) 2 
on U/2D e ff (all energies are to be normalized by the 
effective bandwidth 2 D e ff) 1 which will be confirmed by 
the results of direct numerical computations in the next 
Section (cf. Fig. |U(a)). 

Coefficient B is determined by “square” diagram with 
four Cooper vertices with q = 0, “dressed” in arbitrary 
way by disorder scattering, which is shown in Fig. mb) 

S 


b =? t Y1 < G(ie„;pi,p2)G(-ie n ;-p2,-p 3 ) 


n P1P2P3P4 


G(ie„;p3,P4)G(-*e n ;-p 4 ,-pi) >, 


( 20 ) 


where < • • • > denotes averaging over disorder, while 
G(i£ n \ Pi,P 2 ) (and other similar expressions) represent 
exact single - particle Green’s functions for the fixed con¬ 
figuration of the random potential. Performing standard 
summation over Matsubara frequencies, we obtain: 


B = 


ds , £ 
—— th—— 
2iri 2 T 


^2 < G r {e\ Pi, p 2 ) 


P1P2P3P4 

A ( 


G (-G—P2,-Ps)G A£;p 3 ,p 4 )G (—e;— p 4 ,-pi) > . 


(21) 


Due to zero value of momentum q = 0 in Cooper vertices 
and the static nature of disorder scattering, we can now 
use certain generalization of the Ward identity (fTTH) to 
get (at T = T C ): 


B = 




e/2 T c \ 
ch;2 w- c ) 


N 0 (e) 


( 22 ) 


Detailed derivation is presented in Appendix [A] In BCS 
approximation, using the conduction band of infinite 
width with constant density of states A'o(O), we im¬ 
mediately obtain from Eq. (IXil) the standard result: 

B = ^JVb(O) da. 

Again, replacing here e — > e — /i, to move the origin 
of energy to the middle of the conduction band, we can 
write: 


dE 

J-00 4 (£- m ) 3 



£ — fj, 

~2T7 


(e-n)/2T c \ 

ch^ ) 


iV 0 (£) 


(23) 

It is seen, that disorder dependence of the coefficient 
B (similarly to A) is also determined only by disor¬ 
der widened density of states Nq(e) and chemical po¬ 
tential, so that in the case of semi - elliptic “bare” con¬ 
duction band it is reduced to simple replacement D —> 
D e ff , leading to universal dependence of B(2D e ff) 3 on 
U/2D e ff, which is confirmed by the results of direct nu¬ 
merical computations presented in the next Section and 
shown in Figj4)3. 

It should be stressed, that Eqs. (1X71) and (1XH) for GL - 
coefficients A and B were obtained with the use of exact 
Ward identities, and are thus valid also in the limit of 
strong disorder (beyond Anderson localization). 

Universal dependence on disorder, related to conduc¬ 
tion band widening by disorder scattering, is also valid 
for specific heat discontinuity at T c , as it is completely 
determined by coefficients a and B: 


C s (T c ) 


Cn{T c ) = Tc— 


(24) 


Appropriate numerical results are also given in the next 
Section (cf. Fig. Ob)). 

Coefficient C before the gradient term of GL expan¬ 
sion is determined essentially by two - particle character¬ 
istics (due in particular to non - trivial q - dependence of 
the vertex, which is obviously changed by disorder scat¬ 
tering). In particular, the behavior of G is significantly 
changed at Anderson transition [271] , so that no univer¬ 
sality of disorder dependence is expected in this case. 


IV. MAIN RESULTS 

Let us now discuss the main results of our numerical 
calculations, directly demonstrating the universal depen¬ 
dencies of GL - coefficients A and B and specific heat 
discontinuity at T c on disorder. 
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FIG. 2: Universal dependence of superconducting critical 
temperature on disorder for different values of Hubbard at¬ 
traction. 


In Fig. [2] we show the universal dependence of criti¬ 
cal temperature T c on Hubbard attraction U for differ¬ 
ent levels of disorder, which was obtained and discussed 
in detail in Refs. [20 . [21] ]. Typical maximum of T c at 
U/2D e ff ~ 1 is characteristic of BCS - BEC crossover 
region. 

In Fig. [3] we present disorder dependencies of GL 
coefficients a (Fig. [3[a)) and B (Fig. 03(b)) for different 
values of Hubbard attraction. We can see that a in gen¬ 
eral increases with the growth of disorder. Only in the 
limit of strong enough coupling U/2D >1.4 (curves 4 and 
5) in the region of weak disorder we observe weak sup¬ 
pression of a by disorder scattering. Coefficient B pretty 
fast grows with disorder in the region of weak coupling 
(curve 1 in Fig. 03(b)), while in the region of strong cou¬ 
pling this growth becomes more moderate (curves 4,5 in 
Fig. EKb), so that in this region the dependence of B on 
disorder becomes almost independent of the value of U 
(curves 4 and 5 practically coincide). 

However, this rather complicated dependence of coef¬ 
ficients a and B on disorder is determined solely by the 
growth of effective conduction bandwidth with disorder¬ 
ing given by Eq. ([5]). In Fig. 2] we show the universal 
dependencies of GL - coefficients a (a) and B (b), nor¬ 
malized by appropriate powers of effective bandwidth, 
on the strength of Hubbard attraction. In the absence of 
disorder (dashed line with squares) coefficients a and B 
drop fast with the growth of U. Other symbols in Fig. 
E3 show the results of our calculations for different levels 
of disorder. It is clearly seen, that all the data ideally fit 
the universal curve, obtained in the absence of disorder. 

Coefficients a and B determine specific heat discon¬ 
tinuity at the critical temperature (1M1) . As these coef¬ 
ficients and T c [20l [21 ] depend on disorder in universal 
way due only to the growth of the effective bandwidth 
©, the same type of universal dependence is also valid 
for specific heat discontinuity. In Fig. [5](a) we show the 



FIG. 3: Disorder dependence of GL - coefficients a (a) and 
B (b), normalized by their values in the absence of disorder, 
for different values of Hubbard attraction. 


dependence of specific heat discontinuity dC = C s — C n 
on disorder for different values of Hubbard attraction U. 
It is seen, that in the region of weak coupling (curve 
1) specific heat discontinuity is suppressed by disorder¬ 
ing, for intermediate couplings (curves 2,3) weak disorder 
leads to the growth of specific heat discontinuity, while 
the further growth of disorder suppresses this disconti¬ 
nuity. In the region of strong coupling (curves 4,5) the 
growth of disorder leads to significant growth of specific 
heat discontinuity, which is mainly related to the similar 
growth of T c (cf. 0, Si]). However, this complicated 
dependence of specific heat discontinuity on disorder is 
again completely determined by the growth of the effec¬ 
tive bandwidth ]5]). In Fig]5](b) we show the universal 
dependence of specific heat discontinuity on U, normal¬ 
ized by the bandwidth 2D e ff. Black squares represent 
data in the case of absence of disorder. Other symbols 
in Fig. ]5](b) show the data for different disorder levels. 
It is seen again, that all the data precisely fit the uni¬ 
versal dependence of specific heat discontinuity obtained 
in the absence of disorder. Specific heat discontinuity 
grows with the growth of U in the region of weak coupling 
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U/2D eff 



U/2D e ff 


FIG. 4: Universal dependence of GL - coefficients a (a) and 
B (b) on Hubbard attraction for different values of disorder. 

U/2D e ff <C 1 and drops with the growth of U in the limit 
of strong coupling U/2D e ff 1. The maximum of spe¬ 
cific heat discontinuity is observed at U/2D e ff « 0.55. 
Actually, this dependence of specific heat discontinuity 
qualitatively resembles the similar dependence of criti¬ 
cal temperature, though the its maximum is reached at 
smaller values of Hubbard attraction. 


V. CONCLUSION 

Using the combination of Nozieres - Schmitt-Rink ap¬ 
proximation with the generalized DMFT+E approach we 
have studied disorder influence upon coefficients A and 
B determining the homogeneous Ginzburg Landau ex¬ 
pansion and specific heat discontinuity at superconduct¬ 
ing transition in attractive Hubbard model. 

We have demonstrated analytically, that in the case of 
the “bare” conduction band with semi - elliptic density of 
states disorder influence on GL - coefficients A and B and 
specific heat discontinuity is universal and is controlled 
only by the general conduction band (density of states) 



A/2D 



U/2D e ff 


FIG. 5: Dependence of specific heat discontinuity at critical 
temperature dC = C s — C n on disorder for different values of 
Hubbard attraction U (a) and universal dependence of this 
discontinuity on U for different values of disorder (b). 


widening by disorder scattering and illustrated this con¬ 
clusion with explicit numerical calculations, performed 
for the wide range of attractive potentials U, from weak 
coupling region where U/2D e ff <C 1 and superconduct¬ 
ing instability is described by the usual BCS approach, 
up to the strong coupling region where U/2D e ff 1 
and superconducting transition is determined by Bose - 
Einstein condensation of preformed Cooper pairs. 

These results essentially prove the validity of the gener¬ 
alized Anderson theorem in BCS - BEC crossover region 
and in the limit of strong coupling not only for supercon¬ 
ducting T c [2f| HHj, but also for homogeneous Ginzburg 
- Landau expansion, determining appropriate thermody¬ 
namic effects, like specific heat discontinuity at transition 
point. 

This work is supported by RSF grant 14-12-00502. 
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of momenta and frequencies (i.e. to the replacement 
Gi -O- Gi) in all Green’s functions standing to the right 
from the vertex. Let us assume, that the central of three 
Cooper vertices was inserted in the fourth Green’s func¬ 
tion, as shown in the right part of Fig. (Htb). Arbitrary 
insertion of the first Cooper vertex into one of the first 
four of Green’s functions leads to the following result: 


L \ 


X /■ 


•l \ 


Gi G 2 G 3 G 4 G 5 G s G 7 G b G g 


/ \ 


Ci. L •' 


(b) 


G! G 2 G 3 Gg Gg 


G1G2G3G4 G1G1G2G3G4 + G1G2G2G3G4 + 

G 1 G 2 G 3 G 3 G 4 + G 1 G 2 G 3 G 4 G 4 , (A3) 

so that taking into account G ” 1 — G ” 1 = 2s, we get: 


FIG. 6: Diagrams for coefficient B and derivation of the gen¬ 
eralized Ward identity. 


Appendix A: Coefficient B in the presence of 
disorder 

Coefficient B is determined by “square” diagram with 
four Cooper vertices with q = 0, “dressed” by disorder 
scattering, shown in Fig. mb). Corresponding analytic 
expression was given above in Eq. EDI) . After the stan¬ 
dard summation over Matsubara frequencies B is written 
as in (EH, i.e. is determined by the following combination 
of four Green’s functions with real frequencies: 

< GK ( £ ; Pi; P 2 )G A (s; -P 2 , -Ps) 

P 1 P 2 P 3 P 4 

G R {s; p 3 , p 4 )G A (-e; -p 4 , -pi) > . (Al) 

where < • • • > denotes averaging over disorder and 
G r ( a \s; pi, p 2 ) are the exact retarded (advanced) single 
- particle Green’s functions for the fixed configuration of 
disorder. 

Typical diagram for the fourth order of disorder scat¬ 
tering (dashed lines) is shown in Fig. [DK a )- Arbitrary 
diagrams for such four - particle Green’s function can 
be obtained from diagrams for single - particle Green’s 
function of the same order of disorder scattering by arbi¬ 
trary inserting three Cooper vertices into “bare” electron 
Green’s functions, as shown in Fig. Eta)- Taking into 
account the static nature of disorder scattering and zero 
value of transferred momentum q = 0 in Cooper vertices, 
we can evaluate (ED using certain generalization of exact 
Ward identity (flDl) . derived in Ref. [19,]. 

Let us take the diagram for single - particle Green’ 
function, shown in the left part of Fig. Elb), and consider 
certain configuration of momenta transferred by dashed 
lines. Here we have nine “bare” electron Green’s functions 
with momenta p 4 ■ ■ • pg. In the following we use short 
notations: 

Gi = G R (s; Pi ) Gi = G A (-s; - Pi ), (A2) 

where Gq^ A \s; p) = is the “bare” Green’s func¬ 

tion. Insertion of Cooper vertex leads to the sign change 


G1G1G2G3G4 G1 G i +, 

2s 

s' — 1 sy — 1 

1 '“ r 4 ^4 

~r Or 1 Lt 2 ^3 ^4 ^4 " 


2s 

G P^i S^i XX S~1 S^i XX S~1 

l(jr2Cjr3(_r4 — (jr 10^2^3 *-*4 

2s 


(A4) 


Then G 4 G 5 Gq G 7 G 3 Gg —>■ G 4 G 5 G% G - Gg Gg and after 
all insertions of the last (third) Cooper vertex in one 
of the six Green’s functions G 4 • • • Gg, we again obtain: 

G4G s G fi G 7 G H G fl -G4G s GfiG7G s G 3 

2 e 

Thus we get: 


< G R (s)G A (-s)G R (s)G A (-s) >= 
G a (—s) - G R (s ) G A (-e) - G R (s) 


2s 


2s 


= ^(< G a (-s)G a (-s) > + < G r (s)G r (s ) > - 
—2 < G r (s)G a (-s ) >) = 


-(< G a (-s) >-< G r (s ) >)- 


1 

4 s 2 \ ds 
< G a (—s) >-< G R (s) > 


(A5) 


where we can evaluate two - particle Green’s functions 
with q = 0 again using the analogue of the Ward identity 
(flDl) for real frequencies. Using (IA5I) in (Ell) and making 
in terms with < G A (—s ) > under the integral over s the 
replacement s —> —s, we obtain the following expression 
for coefficient B: 



which was used in the main part of the paper. 
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